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ABSTRACT 

We present for the first time a detailed study of the properties of molecular gas 
in metal-rich environments such as early- type galaxies (ETGs). We have explored 
Photon-Dominated Region (PDR) chemistry for a wide range of physical conditions 
likely to be appropriate for these sources. We derive fractional abundances of the 20 
most chemically reactive species as a function of the metallicity, as a function of the 
optical depth and for various volume number gas densities, Far-Ultra Violet (FUV) 
radiation fields and cosmic ray ionisation rates. We also investigate the response of the 
chemistry to the changes in a— element enhancement as seen in ETGs. We find that 
the fractional abundances of CS, H 2 S, H 2 CS, H 2 0, H 3 0+, HCO+ and H 2 CN seem 
invariant to an increase of metallicity whereas C + , CO, C 2 H, CN, HCN, HNC and 
OCS appear to be the species most sensitive to this change. The most sensitive species 
to the change in the fractional abundance of a— elements are C + , C, CN, HCN, HNC, 
SO, S0 2 , H 2 and CS. Finally, we provide line brightness ratios for the most abun- 
dant species, especially in the range observable with ALMA. Discussion of favorable 
line ratios to use for the estimation of super-solar metallicities and a-elements are also 
provided. 

Key words: astrochemistry - ISM: abundances - ISM:molecules - methods:numerical 
- galaxy:active - stardormation 



1 INTRODUCTION 

Understanding the production of metals, and the subsequent 
enrichment of the ISM is crucial in many field of extra- 
galactic astrophysics. For example, the initial mass func- 
tion (IMF) of stars must have been different in the early 
universe, where no metals were present, and whether there 
is any remaining varia tion of IMF with metallicity remains 
an open question (e.g. lOmukai et al.ll20 05; Santoro fc Shulll 
12009 ; iTumlinsonl |2007l : Ivan Dokkum fc Conrovl l201ll ). Gas 
metallicity is regulated by a complex interplay between star 
formation, infall of metal-poor gas and outflow of enriched 
material. The d iscovery of a r elation between galaxy mass 
and metallicity (|Lequeuxi 1979 h showing that the more mas- 
sive a galaxy is the higher is its metallicity, has increased de- 
bate about the metallicity and its relationship with galaxy 
characteristics and in particular with the properties of its 
interstellar medium. 

Estimating the metallicity of a galaxy remains to date 
a difficult task with few methods available, and these have 



E-mail: bayet@physics.ox.ac.uk; 



large uncertainties. Stellar metallicities are usually derived 
from absorption indices (such as those of magnesium and 
iron), but these weak features are hard to observe at high 
redshift, and only give information about the metallicity 
of t he system in the pas t, when the stars were formed 
(e.g. iMcDermid et al.|[201ll ). Optically determined metallic- 
ities of the ionised gas provide information on the metal 
content of star-forming complexes, allowing astronomers to 
probe the gas phase metallicity further out into the universe 
but this approach only works if the emission spectrum ob- 
served is powered by star formation processes. Where strong 
AGN activity, old stellar populations or shocks dominate, 
this method becomes unreliable. Furthermore, there can be 
strong degeneracies between low and high metallicity popu- 
lations in some ionised gas metallicity indicators. 

In this paper, we explore the possibility of determin- 
ing metallicity via submillimeter observations where no 
dust attenuation is present, remaining a reliable method 
even if star formation does not dominate the ionisation 
of the gas. In particular, we concentrate on exploring the 
high metallicity range represented by early-type galax- 
ies (hereafter called ETGs). These metal-rich objects (el- 
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lipticals and lenticulars) located at the 'end point' of 
galaxy merger sequences, carry important signatures of 
mass assembly and star formation in the Universe (e.g. 
iKavirai et alj 12009). Surprisingly, a significant fraction of 
them (~ 20%) contain some molecular gas res ervoirs and are 
not as 'red-and-dead' as previously thought dWelch fc Sagd 
20031; ISage fc Welch! 12009; ICombes. Young fc Bureau! 120071; 



Sage. Welch fc Youndl2007l ; iKrips et al.ll2010l ; lYoung et all 
201ll : ICrocker et al.ll20l"lF 7 The properties of this molecular 
gas such as its temperature, volume number gas density, col- 
umn density of gas etc. have been estimated (|Bavet et alj 
2012) showing a priori similarity with the molecular gas 
properties seen in the center of our own Galaxy. However, 
the origin of the gas in ETGs remains to date a mystery (but 
see lDavis et alJfeoilbh . Determining the best molecular gas 
metallicity tracer can help to better constrain this origin. 
Indeed, if the metallicity of the molecular gas is similar to 
the measurements obtained from stars and the ionised gas, 
then the origin of the gas in ETGs is likely to be internal due 
to stellar mass loss. Otherwise, it is argued that an external 
origin due to past merging events and/or accretion is more 
likely. 

Here, we perform a first detailed and comprehensive 
theoretical study of the influence of various parameters re- 
lated to metallicity on the molecular gas properties, all likely 
to be appropriate for ETGs. The Photon-Dominated Region 
(PDR) code and the chemical database that we have used 
for our study are described in Section [2 and the physical pa- 
rameters we have selected are discussed in Section[3] Results 
describing the chemistry and its sensitivity to variations of 
these parameters are presented in Section [4] In Section [5l 
we list the best chemical species or combinations of them we 
found which represent well ETG environments. In Section [6] 
we summarize our conclusions. 



2 MODEL DESCRIPTION 

We make use of the UCLJDR model as desc ribed in 
iBell et all |2005l . 120061 ): iBavet et al.l |2009bl , l201ll) . This is 
a time-dependent Photon-Dominated Region (PDR) model. 
In the present application, the code is run for 10 7 yrs, at 
constant density, by which time chemical steady state is 
reached in all models. In fact, the fractional abundances at 
10 6 yrs do not differ significantly from those at 10 7 yrs. Thus, 
time-dependence above 10 6 yrs is unlikely to play a role in 
the chemistry. The code operates in one spatial dimension 
for an assumed semi-infinite slab geometry, and computes 
self-consistently the chemistry and the temperature as func- 
tions of depth and time within the semi-infinite slab, taking 
account of a wide range of heating and cooling processes 
jBavet et al.ll201ll ). In this work, the code is used to deter- 
mine the chemical and thermal properties at all depths up 
to 20 visual magnitudes. 

Here, we improved the UCL_PDR code by adding the 
calculation of the line brightnesses for 13 CO, C ls O, CS and 
HCN transitions (see results in SubsectU). These line bright- 
nesses (in erg s _1 cm" 2 sr -1 ) are calculated us ing the LVG 
approximation (see e.g. Ivan der Tak et ai]|2007l ) and the col- 
lisional rates with H2 are taken from the Leiden At omic and 
Molecular Database (L AMDA; TSchoier et alll2005l ). The de- 
tailed description of these latest changes in the code and 



their consequences in its overall structure are outside the 
scope of this pap e r but can be found in great detail in 
IBell. Bavet fc Vitil (|2012T ). 

The chemical network used here links 131 species in 
over 1800 gas-phase reactions; only H2 is formed by surface 
chemistry. Freeze-out of species on to grain surfaces is not 
considered. This might have a severe impact on the frac- 
tional abundance predictions at A v > 3 — 8 mag since at 
these optical depths, the temperatures, as seen in Fig. [T] are 
low. In addition, freeze-out time is about 10 4 yrs for most of 
the models developed in the present study and 1000 yrs and 
about 100 yrs for Models 13 and 14 respectively. Therefore, 
freeze-out may be important. However, at those extinctions 
the rotational lines considered in this study are mostly opti- 
cally thick in which case freeze-out would not affect the line 
ratios but still the abundances. We have thus run chemical 
models to estimate the impact of freeze-out on abundances 
and this is estimated to be about 3-8%, depending on the 
extinction and species considered. The UCL_PDR code has 
been validated agai nst all other comm only used PDR codes 
and performs well (|Rollig et al.ll2007l ). 



3 PARAMETER SELECTION 

The model requires the setting of a number of physical and 
chemical parameters. The choices we have made are listed 
in Table [fl [2] and 

Most of the calculations are made for an assumed gas 
number density of 1 x 10 4 hydrogen nuclei (in all forms) 
cm -3 , typical value found for molecular gas in galaxies. The 
A v values range from to 20 mags. For Fig. \T\ (left hand 
side) and Fig. [21 we took five representative depths into the 
one-dimensional slab: A v = 1, 3, 5, 8 and 20 mag. These 
values have been selected intending to represent various gas 
phases from a dense PDR surface to dark cloud conditions. 

We have investigated how metallicity changes influence 
the chemistry by studying the effect of an increase of metal- 
licity from 1 to 3.3 zq by steps of 0.2 zq from 1.1 zq. z = 1 zq 
corresponds to solar value of the initial elemental abundance 
ratios (see Table [2]) while z = 2.5 zq means that the solar 
values of the initial elemental abundance ratios have been 
all multiplied by the same factor (of 2.5 in this example). 
The metallicity range has been chosen such as to agree with 
observational results from the ATLAS 3D survey which show 
that metallicity ranges for local ETGs from 0.043 to 0.5 dex 
i.e. from 0.9 z to 3.2 z |McDermid et al.l201ll ). The step 
of 0.2 z© was chosen so as to have a reasonable sampling. 
Several parameters are assumed to scale linearly with metal- 
licity. These are the dust-to-gas mass ratio, the H2 forma- 
tion rate, and the initial elemental abundances. The values 
of these parameters for solar metallicity are given in Table 
[1] and [2] The elemental initial abundances given in Table 
[5] are derived from the most up-to-date esti mates obtained 
on co mets and solar photosphere studies fsee lAsplund et all 
120091 and references therein) . 

High metallicity environments such as ETGs or Milky 
Way metal-rich stars (typically Pop II) have been found 
to have so me of their elements enhanced as compared to 
other s (e.g. iRvan. Norris fc Besselil Il99ll ; iMcDermid et all 
2011). These enhanced elements are oxygen (O), neon (Ne), 
magnesium (Mg), silicon (Si), sulphur (S), argon (Ar), cal- 
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cium (Ca) and titanium (Ti). They are more generally called 
a— elements since their most abundant isotopes are inte- 
ger multi ples of the mass o f the helium nucleus (the a 
particle). lAller fc Greensteinl (|l96Cll ) were the first to see 
a— elem e nt en hancement in metal-poor stars in our Galaxy. 
iPeletierl (| 1989ft was the first to note that you can not repro- 
duce the stellar popul ations of ETG s with o ut a— element 
enhancement whereas iRanalli et all (|2003h : lOrigl iaetaLl 
worked out what this meant in terms of bursty star 
formation. For a review on a — elemen t s and their enhance- 
ments, we refer the reader to [Thomas] l|l999T ). 

In the current UCL_PDR code, only 5/8 a— elements, 
namely O, Mg, Si, S and Ca are included in the chemical 
network (see Table [2] in bold font). For the rest of the pa- 
per, we thus consider only the influence on the chemistry 
of a change in these 5 a— elements (see Subsection 14. 3[) . To 
be more specific, we stud ied the a— element enhancement 
of ISalaris fc Weiss! l| 19981 ) who provided the most recent 
a— elements enhanced mixture able to reproduce well ob- 
servations of the metal-rich globular clusters in the galactic 
disk (see Table [2] column 2) and thought to be appropriate 
for ETGs. In their study, they enhanced the a— elements 
by various factors. We also studied another scenario of 
a— element enhancement used mostly in metal-rich stellar 
evolution model likely appr opriate for ETGs too, investi- 
gated bv IWeiss et alJ (|2004l ). In their paper, contrarily to 
ISalaris fc Weisd (|l998ft . they applied a systematic +0.4 dex 
factor to all the a— elements. Here, we do n ot use directly the 
values reported in ISalaris fc Weissl (|l998l ^ and IWeiss et alJ 
|2004ft because they each use different standard (solar) val- 
ues than ours but we adapted their approach to our standard 
values (see results in Table [SJ. 

The models require the specification of the mean dust 
grain radius and albedo, which affect the transfer of external 
radiation into the cloud. The values given in Table [1] are 
canonical. The external radiation field impinging on the one- 
dimensional semi-infinite slab is assumed to be 1 Habing in 
most cases, however the influence of stronger radiation field 
is investigated through Models 15 and 16 (see Section [4}. 1 
Habing is the standard value used for the Milky Way. The 
cosmic ray ionisation rate standard values used in our study 
is 5 x 10" 17 s" 1 . 

Finally, the models also require the specification of the 
microturbulence velo city for line- width ca lculations . We use 
the value adopted bv lBavet et al.l (|2009bft i.e. 1.5 km s . 

Our grid contains 20 models. Models to 12 investigate 
the effect of metallicity changes on the chemistry. Models 8, 
13 and 14 allow us to study the influence on the chemistry of 
the change of gas number density. Models 15 and 16 look at 
the influence of an enhanced FUV radiation field, whereas 
Models 17 and 18 focus on the influence of the increase of 
cosmic ray ionisation rate. Finally Models 19 and 20 anal- 
yse changes in the high-metallicity chemistry for different 
scenarios of a— element e nhancement. 

As shown in iMcDermid et al.l (|201lft and 
iKuntschner et all |2010ft . some ETGs may have sub- 
solar metallicity stellar populations in their cores. These 
low metallicity cores are likely cau s ed by the accretion of 
low metallicity gas. iKavirai et all |201lft showed that a 
number of ETGs have acquired cold gas through accretion 
during minor m erger events with dw arf galaxies (at redshift 
smaller than D. lDavis et alJ (|2011aT ) also showed that >50% 



Table 2. Initial abundance ratios u sed in Table [3] The a bbre- 
viations 'SW98' a nd 'WE04' refer to lSalaris fc Weisd l|l998l l and 
IWeiss et al.l l l2004ft , respectively (see Sect, [3ll. The standard initia l 
elemental abundance ratios ('ST') are from Asplund et alJ <|200ST l ■ 
Models using the 'ST' values of the initial elemental abundance 
ratios are the models corresponding to z = 1 zq . Bolt font is used 
to represent the a— elements. 





ST 




SW98 




WE04 




Fe/H 


3.16x10" 


-s 


3.16x10" 


-5 


3.16x10" 


-s 


C/H 


2.69x10" 


-4 


2.69x10" 


-4 


2.69x10" 


-4 


O/H 


4.90x10" 


-4 


1.55x10" 


-3 


1.23x10" 


3 


N/H 


6.76x10" 


■5 


6.76x10" 


-5 


6.76x10" 


-5 


Si/H 


3.24x10" 


■5 


8.13x10" 


-5 


8.13x10" 


-5 


S/H 


1.32x10" 


■5 


2.82x10" 


-5 


3.31x10" 


■5 


Ca/H 


2.19x10" 


■6 


6.92x10" 


-6 


5.50x10" 


■6 


He/H 


8.51x10" 


-2 


8.51x10" 


-2 


8.51x10" 


-2 


Mg/H 


3.98x10" 


-5 


1.00x10" 


-4 


1.00x10" 


-4 


Na/H 


1.74x10" 


■6 


1.74x10" 


-6 


1.74x10" 


■6 


Cl/H 


3.16x10" 


-7 


3.16x10" 


-7 


3.16x10" 


-7 



of gas-rich field ETGs in the local volume have obtained 
this gas from mergers or accretion. Low metallicity gas 
that is accreted will, however, quickly be enriched by SNII 
(within few hundred Myrs) from the newly formed stellar 
population. It is thus likely that almost all gas rich ETGs 
have solar metallicity gas, or above. This assumption is 
reinforced by the fact that in cluster environments, where 
most of the ETGs studied in this paper are, the dwarf 
satellite galaxies are coming in fairly dry, having most of its 
gas already stripped out, and so, a larger fraction of the gas 
is actually coming from stellar mass loss and therefore must 
be of higher metallicity. The che mistry of low metal licity 
gas has been partially explored in lBayet~et al. (2009b]) and 
forthcoming more detailed studies are in preparation. 



4 SENSITIVITY OF CHEMISTRY TO 

VARIATIONS OF THE UCL_PDR MODEL 
PARAMETERS 

In this section, we present the trends the chemistry is show- 
ing with respect to the changes of the parameters described 
in Section O We have selected about 25 molecules for closer 
study, either for their likely detectability or for their chemi- 
cal interest. We have arbitrarily fixed the limit of detectabil- 
ity of a molecule X to be with a relative abundance of 
n(X)/n H = 1 x 10" 12 (where n H =n(H) + 2n(H 2 )), a cri- 
terion roughly satisfied in our own galaxy. The rest of our 
paper thus deals specifically with species H 2 , H, electron, 
C+, C, CO, CS, SO, S0 2 ^ OCS, H 2 CS, H 2 S, C 2 H, OH, 
H 2 0, H 3 0+, HCO+, CH, NH 3 , CN, HCN, HNC, H 2 CO, 
H 2 CN and CH3OH. We additionally exclude methanol from 
this study because this species is mainly form on grains and 
such reactions are not included in the model we have used. 
There is only one (minor) gas-phase formation route that we 
include in our models but since we are not properly mod- 
elling its formation, we do not think we can say anything in 
this study about the variation of methanol abundances with 
parameters such metallicity. 
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Table 1. Standard model parameters (see text in Sect. [3). These parameters are those used in Models 0-20 except if a different value is 
specified in Table [3] or described in the text. 



Gas Density 

Mean photon scattering by grain 
External FUV radiation intensity (= 
Cosmic ray ionization rate (=C©) 
H2 formation rate coefficient (=Rq) 
Dust:gas mass ratio (=d0) 
Metallicity (=zq) 
Grain size 
Grain albedo 
A v maximum 
Microturbulent velocity 



I©) 



3x10 



10 4 cm- J 
0.9 

1 Habing" 
5.0xl0~ 17 s" 1 
- 18 x^exp( 7 ^_) em's" 1 
1/100 
solar values b 
0.1 um 
0.7 
20 mag 
1.5 kms -1 



a : The unit of the standard Interstellar Radiation Field (ISRF) intensity is I© = 1.6x10 3 erg cm 2 s" 1 iHabind jl968h : b : z = 1 = 
zq corresponds to solar values of the initial elemental abundance ratios (see Table [2] first column) while z = 2.5 zq means that the 
solar values of the initial elemental abundance ratios have been all multiplied by the same factor (of 2.5 in this example). 



Table 3. Input parameters of the UCL_PDR models used to perform this study (see Sections [3] and [4j. The abbreviation 'ST' represents 
the standard values of the initial elemental abundance ratios while the abbreviations 'SW98' and 'WE04' correspond to the values of 
initial elemental abundance ratios referenced in Table [2] This table does not present all the input parameters of the UCL_PDR code for 
each model but only lists the parameters set to values different from the standard ones (see also Table [T] for complementary information). 



Model 


Metallicity 


Dust-to-gas 


H2 form. 


Ini. Elem. 


FUV rad. field 


C 


Gas Density 




(z©) 


mass ratio (dQ) 


rate coeff. (Rq) 


Abund. ratios 


(I©) 


(C©) 


(cm -3 ) 





1 


1 


1 


ST 


1 


1 


10 4 


1 


1.1 


1.1 


1.1 


STxl.l 


1 


1 


10 4 


2 


1.3 


1.3 


1.3 


STxl.3 


1 


1 


10 4 


3 


1.5 


1.5 


1.5 


STxl.5 


1 


1 


10 4 


1 


1.7 


1.7 


1.7 


STxl.7 


1 


1 


10 4 


5 


1.9 


1.9 


1.9 


STxl.9 


1 


1 


10 4 


6 


2.1 


2.1 


2.1 


STx2.1 


1 


1 


10 4 


7 


2.3 


2.3 


2.3 


STx2.3 


1 


1 


10 4 


8 


2.5 


2.5 


2.5 


STx2.5 


1 


1 


10 4 


9 


2.7 


2.7 


2.7 


STx2.7 


1 


1 


10 4 


10 


2.9 


2.9 


2.9 


STx2.9 


1 


1 


10 4 


11 


3.1 


3.1 


3.1 


STx3.1 


1 


1 


10 4 


12 


3.3 


3.3 


3.3 


STx3.3 


1 


1 


10 4 


13 


2.5 


2.5 


2.5 


STx2.5 


1 


1 


10 s 


14 


2.5 


2.5 


2.5 


STx2.5 


1 


1 


10 6 


15 


2.5 


2.5 


2.5 


STx2.5 


10 


1 


10 4 


16 


2.5 


2.5 


2.5 


STx2.5 


100 


1 


10 4 


17 


2.5 


2.5 


2.5 


STx2.5 


1 


10 


10 4 


18 


2.5 


2.5 


2.5 


STx2.5 


1 


100 


10 4 


19 


2.5 


2.5 


2.5 


SW98 


1 


1 


10 4 


20 


2.5 


2.5 


2.5 


WE04 


1 


1 


10 4 



4.1 Gas number density, FUV radiation field and 
cosmic ray ionisation rate 

As rec ently published by iKrips et al.l (|20T0t ) ; ICrocker et al.l 
l|201ll ). HCN and HCO + lines have been detected in high- 
metallicity environments such as ETGs. These detections 
show the likely presence of gas denser than that traced 
by CO since to excite the transitions HCO + (J=1-0) and 
HCN(J=l-0), part of the gas density needs to be close to 
the corresponding line critical densities of 3.4x 10 4 cm -3 and 
2.3x 10 5 cm -3 , respectively. We investigated thus the change 
in high-metallicity chemistry when the density increased by 
a factor of 10 (Model 13 in Table 0) or 100 (Model 14) as 
compared to the selected standard value of 10 4 cm -3 (Mod- 
els 0). Not surprisingly, changes in gas density affect high- 



metallicity chemistry in the same way as found for solar- 
met allicity c hemistry. Detailed an alysis of such changes can 
be found in iBavet et alj (|2009bT ). Similarly, the influence 
of the FUV radiation field and the cosmic ray ionisation 
rate on high-metallicity chemistry is the same as for solar- 
metallicity chemistry. Since results have been already pub- 
lished, for the FUV radiation field i nfluence we r efer th e 
reader to Figure 5 and Table 8 from IBavet et al.l (|2009bh . 
For a detailed influence of the cosmic ray ionisation rate 
on the chemistry, we refer the reader to e.g . Papadopoulosl 
<|2010l ); IBavet et al.l (|201lD ; iMeiierink et al.l (|201ll ). 
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4.2 Metallicity 

We have plotted in Fig. [T] the temperature distribution ob- 
tained from the self-consistent thermal balance, as a func- 
tion of the metallicity (left hand side) for the five A v selected 
previously. We have also plotted the temperature distribu- 
tion as a function of A v for the range of studied metallic- 
ities (i.e. from solar to 3.3 z© by step of 0.2 z©, see right 
hand side of Fig. [TJ. The curves on the left hand side of 
Fig. [1] show the same general behaviour for A v > 1 mag 
i.e. decreasing with increasing metallicity, reaching a mini- 
mum of about 6 K for A v = 5 mag and going up again for 
higher extinctions. This incre ase of tempe r ature after A v = 
5 mag, already mentioned in lBavet et al.l |201ll ). is due to 
the fact that at high visual extinctions the C and CO lines 
become optically thick and hence are less able to cool the 
gas whilst the other sources of heating (such as the cosmic 
rays) which dominate the heating at these depths (60% of 
the total heating) remains constant. At A v = 5 mag, the 
cooling is actually dominated by CS (35.7%), followed by 
H 2 (20%) and 13 CO (13%) which show again the impor- 
tance of the n ewly added radiat ive cooling mechanisms (see 
Section [21 and iBavet et al]|201ll '). The temperature curve at 
A v = 1 mag appears insensitive to metallicity changes; in- 
creasing the metal content of the ISM thus does not seem 
to have significant impact on the total heating or cooling 
mechanisms of the UV irradiated part of the ISM. This is 
indeed expected since increased metallicity enhances sim- 
ilarly the PAH and dust photo-electric heating as well as 
the CO cooling. The ratio between heating and cooling pro- 
cesses stays thus roughly the same from a metallicity value 
to another. No additional radiative contribution to the total 
cooling from molecules are seen at A v = 1 mag since at such 
depth, molecules are mainly dissociated by the FUV radi- 
ation field. Generally, the temperatures seen in Fig. [T] stay 
below 20 K whatever the metallicity and the A v , ensuring a 
rich chemistry. 

Increasing the metallicity (see Fig. [2]) has a different 
impact from one species to another. When increasing the 
metallicity, we can disentangle three different behaviours in 
the chemical abundances. The species seen on the left hand 
side of Fig. [2] decrease by a factor of 3 to 10 with increasing 
metallicity whereas those plotted on the upper right hand 
side of Fig. [2] increase by a factor of 3 to 15 with increasing 
metallicity. H 2 , C, electron, CS, H 2 S, H 2 CS, H 2 0, H 3 0+, 
HCO + and H 2 CN do not show significant variations and 
thus can be considered as species insensitive to metallic- 
ity increase from the solar value. Interestingly, one can also 
add to this list of 'constant' fractional abundances those of 
OH and CN for the case A v = 1 mag only (see black solid 
lines in Fig(2j. Outside these three categories is SO with an 
increased fractional abundance at A v = 1, 3 and 20 mag 
but otherwise decreasing; S0 2 with a decreased fractional 
abundance at A v = 5 and 8 mag but otherwise increasing; 
NH3 with a decreasing fractional abundances with increas- 
ing metallicity for A v = 3 and 20 mag but otherwise in- 
creasing for A v = 5 mag and increasing then decreasing for 
A v = 8 mag (see bottom right hand side of Fig. [2]). C + and 
C 2 H at A v = 1 mag have smaller variations than those seen 
for SO, S0 2 and NH 3 but the fact that the C + and C 2 H frac- 
tional abundances are high (i.e. above 10 -10 ) makes these 



more moderate variations as meaningful as those seen for 
SO, S0 2 and NH 3 . 

Complementary results on the variation of the fractional 
abundances of chemical species wit h changes in metallicity 
can be found in lBavet et ail |2009bD , where we showed that, 
CS in particular, is sensitive to changes in metallicity. This 
apparent contradiction with our curre nt results ca n be ex - 
plained by the fact that the models in IBavet et al.l |2009b) 
where for low metallicity whereas in the present study we 
investigate su per-solar metallicitie s appropriate for ETGs. 
In addition, in lBavet et~afl (|2009bl ) the FUV radiation field 
is set three times higher than in our present study, which 
has a dramatic effects on the CS fractional abundance (see 
their Table 8). 

It is on average very difficult to determine the reason 
for those changes in the chemistry caused by the metallicity 
because the effects are non-linear involving changes in both 
total elements and dust-to-gas ratio. However, the first cat- 
egory of species can be understood as actually a result of 
what happen to the species of the second category when the 
metallicity increases. For instance, species such as HCN and 
HNC (category 2) form predominantly from CN, a species 
which decreases its abundance at high metallicity (first cat- 
egory). This decrease is linked to the decrease of H at high 
metallicity, H being involved in one of the main formation 
routes of CN. The decrease of the CN abundance is also rein- 
forced by the destruction route CN + H 2 where H 2 remains 
high whatever the metallicity. Thus, if HCN and HNC abun- 
dances increase with metallicity, this must cause the CN 
abundance to drop, and vice-versa. Similar chemical routes 
exist between CO, C 2 H and H 2 CO. Detailed chemical analy- 
sis following the formation and destruction routes of some of 
the species included in our models has been performed and 
confirms the present interpretation of these categories. The 
behaviour of species whose abundances are independent of 
the change in metallicity can be explained by the fact that 
those species are both formed and destroyed by exchange 
reactions that equally depend on metallicity and the net re- 
sult is that no change occurs due to metallicity changes. For 
instance, H 2 S is destroyed mainly by oxygen which increases 
with the metallicity, but is formed mainly by electrons re- 
combining with HsS + , H 2 S + and H 2 Sj, reactions which also 
increase with metallicity. This results in a simultaneous in- 
crease in both formation and destruction, which leaves the 
abundance of H 2 S unchanged. More generally, an increase in 
metallicity leads to a decrease in H, C + and molecular ions 
via dissociative recombination, and an increase in electrons. 
One notes also that, as a secondary effect of the increase in 
metallicity, the secondary photon process is enhanced due 
to an increase of total hydrogen nuclei. 

4.3 a— element enhancement 

This is the first time that high-metallicity chemistry has 
been studied with respect to the changes in a— elements 
content (see Fig. [3J. As described in Section [3] and seen 
in Table [51 we have used in the present work two mix- 
tu res of a— elemen t s, one corresponding to the approach 
of ISalaris fc Weissl (1 19981 ) with different fact ors for each 
q— ele ment (Models 19 in Table[3J) and one from lWeiss et all 
who increase by a factor of 0.4 dex all the a— elements 
(Models 20 in Table EJ). 
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We can divide molecular species in three categories: 
those insensitive to the changes in a— elements i.e. whose 
fractional abundances do not change by more than a fac- 
tor of 2: CO, H, OH, NH 3 , H 2 S, OCS, H 3 0+, HCO+; the 
species which have their fractional abundances increasing 
by a factor between 3 and more than 3 orders of magni- 
tude when the a— elements are enhanced (see right hand 
side of Fig. [3)) i.e. SO, S0 2 , CS and H 2 0. Finally, the third 
category includes CN, HON, HNC, CH, C 2 H, H 2 CS, C+, 
C and H 2 CO. These species decrease in fractional abun- 
dances with the a— element enhancement, by up to 5 orders 
of magnitude for C 2 H (see left hand side of Fig. [3]). We have 
not shown in Fig. [3] the evolution of the fractional abun- 
dance of H 2 CN with respect to the change in a— elements 
because its fractional abundance drops close or below our 
limit of detectability. Hence this species, even if sensitive to 
the changes in a— elements is not likely to be detectable. 
However a stringent non-detection could be indicative. Sim- 
ilarly, for C 2 H, CH at A v > 8 mag, and H 2 CO, the decrease 
of fractional abundance is so pronounced that they fall under 
the assumed limit of detectability. Otherwise, the response 
of the fractional abundance of the species belonging to the 
second and third categories make them quite sensitive, hence 
potential good tracers of a— element enhancement. For those 
species whose fractional abundance is above 10 -10 whatever 
the optical depth, they are likely detectable. 

Another interesting aspect is the potential difference in 
the chemistry between the two scenarios of a— element en- 
hancement. The species showing the most contrast in their 
fractional abundances are C 2 H, H 2 S, CH and C (increase of 
a factor greater than 2 between Models 19 and 20) and SO 
and S0 2 (decrease by a factor greater than 2). Here, only 
species showing a fractional abundance varying by a factor 
of more than five between the scenarios of a— element en- 
ha ncement (i.e. only C 2 H) may be useful for distinguishing 
the|SaI aris fc Weiss! (|l998h scenario (i.e. various factor s ap- 
plied to the a— elements) from the IWeiss et al.l |2004l ') one 
(i.e. same factor of +0.4 dex applied to all a— elements). 
However the fractional abundance of C 2 H in Models 19 and 
20 is very low, whatever the optical depth which makes this 
species unlikely to be detectable. To summarize, it appears 
thus not possible to use models to disentangle the two sce- 
narios of a— element enhancement. 

5 MOLECULAR TRACERS OF 

HIGH-METALLICITY ENVIRONMENTS 

5.1 Fractional abundances 

In this section, we attempt to identify the best chemical 
tracers of high-metallicity environments and to provide cor- 
responding line brightnesses (in erg s" 1 cm' 2 sr" 1 ) useful for 
future observational programmes. We restrict the number of 
species we analyse to those (i) belonging to the list of the 
most sensitive species seen in Subsection 14.21 and 14.31 i.e. to 
those having a fractional abundance changing by more than 
a factor of 5 with changes in metallicity or in a— elements; 
to those (ii) having a fractional abundance at least 10 times 
higher than our arbitrarily fixed limit of detectability; and 
to those (iii) having at least one transition previously ob- 
served in extragalactic environments. This leads to the de- 



tailed study of C+, C, CO, C 2 H, CN, HCN, HNC, OCS, CS, 
SO, S0 2 and H 2 0. 

From Subsection 14.21 species likely to be the best at 
probing changes in metallicity are C + , C 2 H, CN, HCN, 
HNC, OCS, CO. Similarly, the most likely chemical trac- 
ers for determining whether there is an enhancement in 
a-elements are C+, C, CN, HCN, HNC, SO, S0 2 , H 2 
and CS. Potentially, a non-detection of C 2 H may also be 
used to indicate the presence of an enhanced a— elements 
environment. 

In Table[4] we provide values of fractional abundance ra- 
tios involving the above species with CS and CO for Models 
0,8 and 12 at 5 different extinctions. The fractional abun- 
dance of CS and CO have been selected as denominators for 
the ratios because CS is shown to be rather insensitive to 
changes in metallicity and CO insensitive to the presence of 
an ot — element enhancement (see Section [4|. 

There is no strong variation seen in the ratios for A v ^ 3 
mag. Therefore, this range of extinctions does not have an 
observational interest. The ratios showing the largest vari- 
ations when the metallicity increases are C + /CO, SO/CO, 
CN/CO and CN/CS, with a factor of differences between 
Models 0, 8 and 12, at A v = 5 mag, of 160, 36, 22 and 
16, respectively. These ratios thus might be considered as 
the best and most likely observable probes of the metallicity 
changes in metal-rich environments. In addition, the ratios 
showing the largest variations when there is an a— element 
enhancement are C 2 H/CS, S0 2 /CO, SO/CO, HCN/CS and 
HNC/CS with factors of difference between Models 8 and 
19, at A v = 5 mag, of >5000, 1049, 669, 117 and 71, respec- 
tively. At higher A v , C/CO becomes also a good indicator 
as it differs between Models 8 and 19 by more than a factor 
of 55. These ratios might thus be considered as the most 
reliable probes of an a— element enhancement in metal-rich 
environments. 

5.2 Line brightness ratios 

The LVG treatment used in the UCL_PDR code to calculate 
the line brightness requires collisional rates for the studied 
species (typically with H 2 ). Unfortunately, for C 2 H, these 
collisional rates are not available. We therefore do not have 
line brightness predicted by our models for this species. We 
also do not calculate line brightness for S0 2 and H 2 be- 
cause they contain a very large number of level populations 
which are required for a proper treatment, and, as a conse- 
quence, an unreasonable amount of computational time. 

In Fig. 21 for Models and 12, we have therefore plot- 
ted line brightness ratios for C+/CS, CN/CS, HCN/CS, 
HNC/CS, CO/CS and OCS/CS as a function of A v . We 
restricted our plots to line ratios involving transitions ob- 
servable simultaneously using large backends on telescopes 
such as IRAM-30m or ALMA, i.e. to those with close 
frequencies such as HCN(l-0)/CS(2-l), HNC(3-2)/CS(5-4), 
OCS (28-27) /CS (7-6), etc. In the right hand side of Fig. |1 
we see the distribution of the line brightness ratios with 
respect to A v for Model only. On the left hand side of 
Fig. [4] we see the comparison between Models (black 
lines) and 12 (grey lines) i.e. the comparison between so- 
lar and 3.3 zq metallicity, for only a selected number of line 
brightness ratios. In Fig. [5] we have plotted the C + /CO, 
C/CO, CN/CO, HCN/CO, HNC/CO and CS/CO ratios. 
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Figure 1. Temperature profiles with respect to Left: the metallicity (in zq) and Right: the optical depth A v (in mag). On the right hand 
side plot, vertical black lines refer to the 5 A v selected and seen on the left hand side plot (see Section [3}. We kept the same line style in 
both plots as well as in Fig.[2]ic. A v of 1 mag is represented by solid lines, A v of 3 mag is shown by dotted lines, A v of 5 mag is designated 
by dashed lines, A v of 8 mag is represented by dashed-dotted lines and A v of 20 mag is described by black dotted-dotted-dashed lines. 
The curves seen on the right hand side plot are obtained for various metallicity (from solar to 3.3 solar by steps of 0.2), using a different 
line style for each metallicity value from 1 solar (solid lines) to 3.3. solar (large dashed lines). 



Similarly to figure [4] on the right hand side of Fig. [5] the 
results from Model 8 are displayed whereas on the left hand 
side of Fig. [5] we see the comparison between Model 8 
(black lines, without a— element enhancement) and Model 
19 (grey lines, with a— element enhancement). The top pan- 
els in both figures present the ratios between well-known 
ISM cooling lines i.e. C+, 01 (146 ^m), 01 (63 ^m) , 01 
(44 pm), CI (610 ^m), CI (370 ^m) and CI(230 ^m) and 
the CS(2-1) line (Fig|J) or the CO(l-O) transition (Fig [5]). 
The CS(2-1) and CO(l-O) transitions were selected as de- 
nominators firstly because the CS and CO fractional abun- 
dances are insensitive to changes in metallicity and presence 
of a— elements, respectively, and secondly, because these 
two lines have easily detectable emission at 3mm in exter- 
nal galaxies (see e . g. Ilsrael fc Baasll200ll . 120031 : iBavet et al.l 
120041 . 12006I . l2009al ; lAladro et al.l 12011^ One notes that OI 
is not amongst the species respecting conditions i), ii) and 
iii), but we have nevertheless plotted its OI (146 /jm), 
OI (63 (im) and OI (44 /zm) line brightnesses because 
they are routinely observed in extragalactic environments. 
Additional line ratios will soon be available online (see 
|http://www.homepages.ucLac.uk/^ucapdwi/interface| ). 

The results derived from figures [4] and [5] can be summa- 
rized as follow: 

• The increase in metallicity impacts similarly for all the 
ratios involving the well-known ISM cooling lines with re- 
spect to the CS(2-1) line, i.e. these ratios have higher values 
up to a factor of 20 (for OI (44 /im)/CS(2-l)) as compared 
to their value in the case of solar metallicity (i.e. Model 0). 
Similarly, the OCS/CS line brightness ratios increase with 
the increase of metallicity, whatever the transition consid- 
ered. For well-known ISM cooling line ratios, the increase 
between solar and 3.3 zq is constant with respect to the A v 
whereas for the OCS/CS ratio, this increase decreases when 
A v increases. This allow us to propose this ratio as powerful 
diagnostics of either the metallicity or the opacity. Indeed, 
if the OCS(8-7)/CS(2-l) line brightness ratio is measured in 
external galaxies to be for instance 0.1, this corresponds ei- 
ther to a gas metallicity of 3.3 solar and a rather translucent 
gas (with an A v ~ 1.8 mag) or the metallicity of the gas is 



solar and the A v ~ 4 mag. To break this metallicity-opacity 
degeneracy, independent estimates of the A v are required. 
To do so, molecular isotopologues detections are essential. 
Possible way to break this degeneracy is to use a set of line 
ratios, providing they involve lines with similar critical den- 
sities, allowing us to assume they originate from the same 
gas component. For example, a combination of H2CS(7-6), 
HCO+(3-2), H 2 CO(3-2), SO(5-4) and C 2 H(3-2) compared 
with CS(5-4) would be useful since those lines, close in fre- 
quency to the CS(5-4) transition, have also similar critical 
densities (of about 3 — 8 x 10 5 cm -3 ). For the other ra- 
tios seen in Fig. H (i.e. CN/CS, HCN/CS, HNC/CS and 
CO/CS), but excluding OCS(8-7)/CS(2-l) the metallicity 
boost to 3.3 z© increases the ratios as compared to the solar 
case for A v ^3 — 4 mag up to a factor of 8 for HCN(1- 
0)/CS(2-l) and HCN(4-3)/CS(7-6). Then, at A v ~ 5 mag 
(A v ~ 2.8 mag for the CN_1/CS(2-1)), the two metallicity 
models give similar ratios. At these extinctions, it will thus 
not be possible to use the models predictions for estimating 
the gas metallicity. Finally, at A v ^ 5 mag, the ratios for the 
solar and 3.3 zq metallicities are decoupled again, with lower 
ratio values for the 3.3 zq metallicity case than for the so- 
lar one. At high extinctions, however the difference between 
the two modelled ratios might not be large enough to be de- 
tectable, leading to a rather uncertain estimate of the metal- 
licity. Here, we have considered observational uncertainties 
(calibration, pointing, etc) of 20%, i.e. to be detectable, we 
need a difference between ratios greater than a factor of 5. 
Two exceptions stand in the HCN/CS and HNC/CS ratios 
involving transitions at about 300 GHz frequency: if the A v 
is estimated independently to be of about 8 mag, then the 
metallicity can be determined (factor of 8 difference between 
the ratios at solar and 3.3 solar metallicity). One notes that 
line brightness CO/CS ratios are not the best ratios to use 
for estimating the metallicity since the difference between 
the two model predictions are not as large as those seen for 
the HNC/CS and HCN/CS ratios (the latter ones being a 
much better indicator of metallicity). 

• In high-metallicity environments such as ETGs (here we 
assume the metallicity to be at 2.5 zq - see Fig- [5j , it may 
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Figure 2. Fractional abundances (n(X)/nn) of species X the most sensitive to the changes in metallicity, in logarithmic scale, where 
njj is the total number of hydrogen atoms, with respect to the metallicity (in z®, see text in Subsect. I4.21 ). These plots summarize the 
chemistry changes seen for Models to 12 with A v of 1 mag, 3, 5, 8 and 20 mags (see caption of Fig. [TJ- The lines for A v of 1 mag in 
the NH3 and the SO2 do not appear because they are well below our limit of detectability. On the left and right (top) hand sides of this 
figure, the species fractional abundances decrease and increase, respectively, with the metallicity changes. On the bottom right hand side 
are the outsidders species behaving differently from an A v to another in response to the metallicity changes. 



possible to confirm the presence of a— element enhancements 
by measuring line brightness SO/CO ratios. Indeed these ra- 
tios show the strongest increase of values between Models 8 
and 19, especially for A v ^ 4 mag and this, whatever the 
transition studied. An a— element enhancement decreases 
the rest of the line brightness ratios by a factor 0.5 to 5. One 



notes that the decrease of the ratios in the case of low fre- 
quency ratios such as HCN(1-0)/CO(1-0), HNC(1-0)/CO(1- 
0) and CN_1/CO(1-0) is not large enough for ensuring a con- 
firmation of a— element enhancement whereas for the ratio 
involving higher frequency line brightnesses such as HCN(4- 
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3)CO(3-2), HCN(4-3)/CO(3-2), especially for A v > 5 mag, 
this is becoming possible. 

In Figure [6] we plot theoretical and observed CO line 
ratios for ETGs. CO(l-O), CO(2-l), 13 CO(1-0), 13 CO(2- 
1), HCO+(1-0) and HCN(l-O) data ha ve been published 
for a limited number of gas- rich ETGs l|Krips et al.l l2010l : 
lYoung et alll201lf ). From the 12 CO/ 13 CO ratio and assum- 
ing a Milky Way value of the abundance ratio 12 C/ 13 C=70 
(see iMartfn et al.l l2010j and references therein), using the 
Eq. 8 in Bavet et al.l l|2004T ). we have obtained an A v for 
the 12 CO gas of 5.5 mag, 8.0 mag, 22 mag and 23.7 mag for 
NGC4150, NGC3032, NGC4526 and NGC4459, respectively 
from 12 CO(1-0)/ 13 CO( 1-0) ratios of 13.6 7, 9. 46, 3.46 and 
3.20, respectively from iKrips et al.l |2010| and lYoung et al.l 
120111 . 

On the left hand side of Fig. |6j we compare the CO 
beam-dilution corrected observations to our model pre- 
dictions. In this plot, we have plotted the variations of 
the CO(2-1)/CO(1-0) line brightness ratio with respect 
to A v for a solar (black line) and a 3.3z© metallicity 
(grey lines). Crosses represent the current observed CO (2- 
1) /CO (1-0) line brightness ratios for four observed gas-rich 
ETGs (derived from the beam c orrected velocity-i ntegrated 
line intensity ratio s displayed in lKrips et al.ll2~010l data. For 
lYoung et al.l l201ll d ata, we used beam dilution correction 
method described in ICrocker et aTll201lD . We can see that 
the variation of the CO(2-1)/CO(1-0) ratio with the increase 
of metallicity is not significant enough to reliably constrain 
the metallicity of these gas-rich ETGs, even if the A v has 
been determined using 13 CO (1-0) data. However, one sees 
that the line brightness CO(6-5)/CO(l-0) ratio can be use- 
ful as a CO line ratio alternative to estimate metallicity. The 
typical erro r on each observed CO(2-1)/CO(1-0) ratio is of 
~ 12% fsee lKrips et~aLll2010l ; lYoung et al.ll201lli . If similar 
error can be obtained on the line brightness CO (6-5) /CO (1- 
0) ratio, then, from Fig|6j reliable metallicity estimate may 
be possible to obtain, especially if the gas is strongly opaque 
as in NGC4526 and NGC4459 (whose Av is outside the 
plot). CO(J-J-1)/CO(1-0) lines at high frequency seem thus 
a better indicator of the metallicity in high metallicity en- 
vironments than lower frequency CO line ratios, but, in any 
case, ratios such as HCN/CS and HNC/CS are more power- 
ful diagnostics. Future work (|Davis et al.ll2012T ) will explore 
this issue with new a nd recently released observational data 
l|Crocker et al.ll201ll ). 



iBavet et all (l2009bh ; iPapadopoulol i|2010l ); iMeiierink et all 
( 201ll ); lBavet et al.l(|201ll ). 



From our study, it appears that the best observable 
tracers of metallicity are C+, C 2 H, CN, HCN, HNC, OCS 
and CO whereas the most likely observable chemical trac- 
ers for probing enhancement in a— elements are C + , C, 
CN, HCN, HNC, SO, S0 2 , H a O and CS. Fractional abun- 
dances and line brightness ratios with respect to CS (in- 
sensitive to metallicity changes) and to CO (insensitive to 
a— element changes) provide more quantitative values, use- 
ful for future follow-up observational programmes. Power- 
ful tracers of metallicity changes have been identified (from 
the most sensitive to the least): the HCN/CS and HNC/CS 
line brightness ratios at all frequency for a gas with A v ei- 
ther smaller than 3 mag or of about 8 mag, line brightness 
OI/CS(2-l) ratios for a gas with A v < 4 mag and OCS/CS 
line brightness ratios at all frequencies for gas with A v ^ 4 
mag. The ideal line brightness ratio for probing the enhance- 
ment of a— elements is the SO/CO whatever the optical ex- 
tinction and the frequency of the lines involved. One notes 
that ratios involving two CO lines are not the best probes 
of either metallicity or a— element enhancement. However 
when CO(6-5)/CO(l-0) line brightness ratio is used and the 
error on the ratio is small, some reasonable estimate of the 
metallicity can be obtained. This however may be difficult 
without a good treatment of beam dilution effects. 

To conclude, if CS data are obtained in sources like 
ETGs, HCN/CS line brightness ratios can be calculated and 
from the results displayed here, we may be able to estimate 
the molecular gas metallicity in these sources which could 
put strong constraints on the origin of the gas in ETGs. This 
work can also be useful in predicting metallicity indepen- 
dently from previously used (i.e. mainly optical) methods in 
other high metallicity environments such as galaxy centres, 
etc. Telescopes such as ALMA and the IRAM-30m telescope 
are amongst the most privileged instruments to use for this 
study since large backends are available, hence observing 
simultaneously several molecular lines of CS, HCN, etc is 
possible. 



6 CONCLUSIONS 

The main results of this theoretical study are contained in 
Section [3] which shows how the chemistry evolves with re- 
spect to the changes in various parameters likely to be ap- 
propriate for high-metallicity environments such as ETGs: 
FUV radiation field, density, cosmic-ray ionisation rate, 
metallicity (increase up to 3 times solar), and an enhance- 
ment in the a— elements. The impact of the last two param- 
eters on the chemistry have never been investigated so far 
in such detail. Table [4] and Figs. 12161 are thus of particular 
interest, especially when considering the fact that at high 
metallicity the influence of the FUV radiation field, density 
and cosmic-ray ionisation rate on the chemistry does not dif- 
fer significantly from what has been previously shown by e.g. 



Table 4. Fractional abundance ratios which are likely the best tracers of high-metallicity environments for various models (see characteristics in Tables [T] to (3J . 



ratio 


C+/CS 


C 2 H/CS 


CN/CS 


HCN/CS 


HNC/CS 


OCS/CS 


co/cs 


C+/CO 


C/CO 


CN/CO 


HCN/CO 


HNC/CO 


SO/CO 


S0 2 /CO 


H 2 0/CO 


A v — 1 mag 
Model 
Model 8 
Model 12 
Model 19 


1.40(1) 
2.87 
3.15 
2.76 


4.97(-l) 
2.10(-1) 
3.00(-l) 
1.00(-2) 


3.98 
1.31 
1.26 
6.50(-l) 


5.57 
9.69 
1.11(1) 
1.29 


3.44 
4.49 
4.81 
5.50(-l) 


4.10(-1) 
1.53 
1.37 

7.30(-l) 


2.68(5) 
3.50(5) 
2.85(5) 
2.32(5) 


5.23(-5) 
8.21(-6) 
l.ll(-5) 
1.19(-5) 


3.70(-l) 
2.50(-l) 
4.20(-l) 
2.30(-l) 


1.49(-5) 
3.75(-6) 
4.42(-6) 
2.81(-6) 


2.08(-5) 
2.77(-5) 
3.91(-5) 
5.54(-6) 


1.29(-5) 
1.28(-5) 
1.69(-5) 
2.39(-6) 


5.79(-7) 
2.58(-7) 
1.96(-7) 
1.32(-6) 


9.76(-10) 
1.20(-9) 
1.07(-8) 
9.98(-9) 


1.41(-5) 
6.47(-6) 
5.17(-6) 
3.05(-5) 


A v = 3 mag 
Model 
Model 8 
Model 12 
Model 19 


2.1(-2) 
5.07(-3) 
4.71(-3) 
1.37(-3) 


1.02(-2) 
2.06(-4) 
1.48(-4) 
1.33(-6) 


4.26(-2) 
9.33(-3) 
7.50(-4) 
6.20(-3) 


3.14(-1) 
7.19(-1) 
8.36(-l) 
1.16(-2) 


2.06(-l) 
3.29(-l) 
3.52(-l) 
5.23(-3) 


2.06(-l) 
1.05 
1.61 

1.22(-1) 


2.81(3) 
4.98(3) 
6.96(3) 
7.55(2) 


7.40(-6) 
1.15(-6) 
6.76(-7) 
1.81(-6) 


1.23(-1) 
5.23(-2) 
4.19(-2) 
3.62(-2) 


1.52(-5) 
1.87(-6) 
1.08(-6) 
8.21(-6) 


1.12(-4) 
1.44(-4) 
1.20(-4) 
1.54(-5) 


7.34(-5) 
6.61(-5) 
5.06(-5) 
6.92(-6) 


4.29(-6) 
1.67(-6) 
1.29(-6) 
5.37(-5) 


4.30(-7) 
4.56(-7) 
3.77(-7) 
2.91(-5) 


9.56(-4) 
5.42(-4) 
3.55(-4) 
2.99(-3) 


A v — 5 mag 
Model 
Model 8 
Model 12 
Model 19 


1.32(-3) 
3.54(-4) 
2.15(-4) 
6.27(-5) 


5.83(-5) 
1.34(-5) 
7.23(-6) 
2.53(-9) 


4.10(-3) 
5.62(-4) 
2.62(-4) 
6.60(-4) 


4.52(-2) 
1.17(-1) 
1.19(-1) 
1.00(-3) 


3.17(-2) 
5.71(-2) 
5.37(-2) 
8.01(-4) 


8.68(-2) 
5.07(-l) 
7.21(-1) 
4.27(-2) 


2.25(2) 
3.26(2) 
3.18(2) 
5.39(1) 


5.89(-6) 
1.09(-6) 
6.78(-7) 
1.16(-6) 


9.54(-2) 
5.27(-2) 
4.61(-2) 
5.46(-3) 


1.82(-5) 
1.72 (-6) 
8.23(-7) 
1.23(-5) 


2.01(-4) 
3.58(-4) 
3.72(-4) 
1.86(-5) 


1.41(-4) 
1.75 (-4) 
1.69(-4) 
1.49(-5) 


8.38(-6) 
9.56(-7) 
2.35(-7) 
6.40(-4) 


4.55(-6) 
2.03(-6) 
7.73(-7) 
2.13(-3) 


7.31(-3) 
5.19(-3) 
3.88(-3) 
2.94(-2) 


A v = 8 mag 
Model 
Model 8 
Model 12 
Model 19 


7.25(-4) 
2.58(-4) 
1.71(-4) 
3.59(-5) 


2.50(-5) 
8.83(-6) 
5.80(-6) 
5.03(-ll) 


2.83(-3) 
5.92(-4) 
2.71(-4) 
7.78(-5) 


2.23(-2) 
7.57(-2) 
9.20(-2) 
2.62(-4) 


1.56(-2) 
3.64(-2) 
4.13(-2) 
1.93(-4) 


4.75(-2) 
3.41(-1) 
5.50(-l) 
1.63(-2) 


1.39(2) 
2.56(2) 
2.61(2) 
3.48(1) 


5.19(-6) 
1.01(-6) 
6.53(-7) 
1.03(-6) 


7.28(-2) 
4.61(-2) 
4.33(-2) 
7.75(-4) 


2.03(-5) 
2.31(-6) 
1.03(-6) 
2.23(-6) 


1.60(-4) 
2.96(-4) 
3.51(-4) 
7.52(-6) 


1.12(-4) 
1.42(-4) 
1.58(-4) 
5.53(-6) 


1.58(-5) 
2.60(-6) 
7.71(-7) 
6.02(-3) 


7.68(-6) 
4.68(-6) 
2.24(-6) 
1.68(-2) 


8.66(-3) 
5.48(-3) 
4.68(-3) 
3.61(-2) 


A v — 20 mag 
Model 
Model 8 
Model 12 
Model 19 


6.74(-4) 
1.65(-4) 
1.16(-4) 
3.37(-5) 


2.05(-5) 
4.10(-6) 
2.97(-6) 
3.72(-ll) 


2.63(-3) 
6.02(-4) 
3.68(-4) 
4.15(-5) 


2.11(-2) 
3.80(-2) 
4.83(-2) 
1.27(-4) 


1.47(-2) 
1.76(-2) 
2.07(-2) 
4.91(-5) 


3.61(-2) 
1.82(-1) 
3.03(-l) 
l.U(-2) 


1.38(2) 
1.91(2) 
2.12(2) 
3.27(1) 


4.88(-6) 
8.64(-7) 
5.49(-7) 
1.03 (-6) 


6.20(-2) 
3.38(-2) 
3.17(-2) 
6.54(-4) 


1.90(-5) 
3.14(-6) 
1.74(-6) 
1.27(-6) 


1.55(-4) 
1.99(-4) 
2.27(-4) 
3.88(-6) 


1.06(-4) 
9.21(-5) 
9.77(-5) 
1.50(-6) 


1.89(-5) 
7.12(-6) 
3.93(-6) 
7.59(-3) 


7.20(-6) 
9.41(-6) 
7.84(-6) 
1.63(-3) 


8.53(-3) 
5.42(-3) 
4.63(-3) 
3.57(-2) 
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Figure 4. Line brightness ratios for various species with respect to the CS lines (study of the metallicity) . We present plots from top 
to bottom of the far-infrared fine structure cooling lines : C+/CS(2-1), OI(44/zm)/CS(2-l), OI(63Atm)/CS(2-l), OI(146/im)/CS(2-l), 
CI(230/W)/CS(2-1), CI(370/im)/CS(2-l), CI(610/im)/CS(2-l), CN/CS, HCN/CS, HNC/CS, CO/CS and OCS/CS ratios. The CNi, 
CN2 and CN3 transitions represent emission at 113.490 GHz, 226.874 GHz and 340.247 GHz, respectively. Every line involved in a ratio 
is mentioned on the right and left hand side plots. Left: Comparison between Models (black lines) and Model 12 (grey lines). Right: 
Predictions for Model 0. The ratio plotted on the left hand side are the ones plotted, in the same line style, on the right hand side. 
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Figure 5. Line brightness ratios for various species with respect to the CO lines (study of the a— element enhancement). We present 
plots from top to bottom of the far-infrared fine structure cooling lines : C+/CO(1-0), OI(63/im)/CO(l-0), OI(146/im)/CO(l-0), 
CI(370/*m)/CO(l-0), CI(610/im)/CO(l-0), CN/CO, HCN/CO, HNC/CO, SO/CO and CS/CO ratios. The CNi, CN 2 and CN 3 tran- 
sitions correspond to the CN emission at 113.490 GHz, 226.874 GHz and 340.247 GHz, respectively. The numbers 1 to 7 in the plot 
of the HCN/CS line brightness ratios correspond respectively to HCN(1-0)/CO(1-0), HCN(2-1)/CO(1-0), HCN(3-2)/CO(2-l), HCN(4- 
3)/CO(3-2), HCN(5-4)/CO(4-3), HCN(6-5)/CO(5-4) and HCN(7-6)/CO(6-5). Similar lines ratios are plotted for HNC/CS. For the plot 
of the SO/CO, we have the SO emission at 1=109.252 GHz/CO(l-0), 2=246.404 GHz/CO(2-l), 3=336.553 GHz/CO(3-2), 4=461.755 
GHz/CO(4-3) and 5=609.96 GHz/CO(6-5). Finally for the CS/CO plot, wc have 1=CS(1-0)/CO(1-0), 2=CS(2-1)/CO(1-0), 3=CS(3- 
2)/CO(l-0), 4=CS(4-3)/CO(2-l), 5=CS(5-4)/CO(2-l) and 6=CS(7-6)/CO(3-2). Left: Comparison between Models 8 (black lines) and 
19 (grey lines). Right: Predictions from Model 8. The ratio plotted on the left hand side are the ones plotted, in the same line style, on 
the right hand side. 
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Figure 6. Example of results obtained when comparing observed and predicted CO line ratios. Right: Results from Model (solar 
metallicity); Left: Comparison between Model (black lines) and Model 19 (3.3 zq - grey lines) predictions for the CO(2-1)/CO(1-0) 
(solid lines) and the CO(6-5)/CO(l-0) (dashcd-dottcd-dot ted lines) line br i ghtness ratios. On top of which we added observational values 
(see black crosses) which represent the data presented in lKrips et alj[2010HYoung et alj|201ll . taking into account of their errors (crosses 
show error bars). The A v derived for NGC4459 (23.7 mag) is outside the plotted range of optical depths hence does not appear on the 
figure. 



